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THE SZEKERES CLASS OF VACUUM 
SOLUTIONS 

The first exact solution to be published which describes a collision be- 
tween plane waves was in fact that of Szekeres (1970). This describes the 
collision of plane gravitational waves with step wavefronts. The solution 
discussed in Chapter 3, which describes the collision between impulsive 
gravitational waves, was published a little later by Khan and Penrose 
(1971). Following this, a substantial paper in which the subject is anal- 
ysed in detail was produced by Szekeres (1972). This paper includes a 
derivation of the field equations in the form given in Chapter 6, a general 
class of exact solutions which includes the two mentioned above as special 
cases, and a discussion of the singularities that arise. The purpose of this 
chapter is to discuss the properties of this general class of solutions. 

9.1 The solution in region IV 

The vacuum field equations appropriate to region IV of the colliding 
plane wave problem may be taken here in the form of equations (6.22a-/). 
These are second order differential equations for the four metric functions 
U(u,v), V(u,v), W(u,v) and M(u,v). 

Equation (6.22a) may immediately be integrated as in (6.24) to give 

e~ u = f + g (9.1) 

where / = f(u) and g = g(v) are arbitrary decreasing functions which, 
according to the appropriate boundary conditions (7.3), are required to 
satisfy /(0) = ±, f (0) = and g(0) = ±, g'(0) = 0. 

In the Szekeres class of solutions, the approaching waves have con- 
stant aligned polarization. In these solutions therefore 

W = 0, (9.2) 

and the main equations (6.22a 1 , e) reduce to the single equation 

2V W = U U V V + U V V U . (9.3) 
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64 The Szekeres class of vacuum solutions 
For this, Szekeres has obtained the solution 

/ 1 - f\ 1/2 /I - \ 1/2 

F = -2/ci tanh -1 yiT^TgJ ~ 2k * tanh -1 ^f^J J ( 9 - 4 ) 

which contains two arbitrary constants 1 k\ and /c 2 - The expression (9.4) 
may also be written in the alternative form 



(9.5) 



With this expression for V, the remaining equations in (6.22) may 
be integrated to give 

M = - log(c/ g ) ^ log(/ + g) 

+ % log(| - /) + | log(| + /) + I log(i - <?) + % log(| + <?) 




+ 2M 2 log yi^f^g + y/l + fy/i+gj 



(9.6) 



where c is a constant. 

As anticipated in the discussion in Section 7.2, it may be observed 
that this expression contains the necessary multiples of log(^ — /) and 
log(| — g) that are required to cancel the effects of the unbounded term 
log/V on the boundary. If the leading terms in the expansions for / and 
g are 

/=|-(c lU ) ni + ..., g=\-( C2 v)^ + ... (9.7) 
where rii > 2, i = 1, 2, then e~ M is continuous across the boundaries if 

h 2 = 2(l-l/m), /c 2 2 = 2(1 - l/n 2 ). (9.8) 

It may thus be observed that the constants k\ and /c 2 are restricted to 
the range satisfying 

1 < k 2 < 2, i = 1,2. (9.9) 



It may be noticed that the constants used here are minus one half of those 
used by Szekeres. 
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It is also appropriate to choose 

c = (cin l c 2 n 2 )~ 1 (9.10) 

to achieve the usual flat metric (3.6) in region I. 

When evaluating the components of the Weyl tensor, it is convenient 
to introduce the new function 



H = 2h yj\ - f y/ \ + 9 - 2k 2 y/ \- gyj \ + f. (9.11) 

With this, it can be seen that 

Vu = -y^- g H f> v v = jTg- 11 ^ ( 9 - 12 ) 

and the non-zero Weyl tensor components are 

*5 = "i(/^(l -H,H„) (9.13) 

In his paper, Szekeres chooses 

/=i-(ciu) ni and g=l-(c 2 v) n > (9.14) 

exactly. In this case, the terms ^ log(| _ /) and ^ log(^ — g) in (9.6) 
exactly cancel the term \og(cf'g'), and M is given by 



-M 



(J _|_ ^(fe?+fei+2fciA:2-l)/2 



2fcifc 2 



(9.15) 

This form is still completely general since it has, in effect, simply used the 
coordinate freedom (6.7). However, it is sometimes convenient to retain 
this freedom, so this restriction will not be made in the following section. 

It may now be pointed out that, with the restriction (9.14), the above 
solution includes the Khan-Penrose (1971) solution for colliding impulsive 
waves when n\ = n 2 = 2, k\ = k 2 = 1. It also includes the Szekeres (1970) 
solution for colliding step waves when n\ = n 2 = 4, k\ = k 2 = a/3/2. 
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The Szekeres class of vacuum solutions 



9.2 The approaching waves 

Having obtained an exact solution in region IV, the question now is to find 
the initial conditions which give rise to it. In fact, it is quite simple to use 
the method described in Section 7.2 to extend any solution in region IV 
back into regions II, III and I. We obtain the corresponding solutions in 
regions II and III simply by replacing g by \ and / by | alternately. In this 
section we will concentrate on region II, and put U = — log sf \ + f. The 
equivalent solution in region III can be obtained from this by replacing / 
by g, and by interchanging u and v. 

Using the suggested method, it can be seen that the solution in region 
II must have the line element 

ds 2 = 2e~ M dudv -(§ + /) (e v dx 2 + e~ v dy 2 ) (9.16) 

where, retaining the coordinate freedom in u, 

/ = l-(ci«r + ... 



1 + \ fV~~f) ' 2 ni (9 - 17) 



f (l +/ )(fcM)/2 

cinxd-/)*?/ 2 ' 



In order to interpret this solution as a plane wave, it is appropriate 
to transform the metric to the Brinkmann form 



ds 2 = Idudr + h{u)(x 2 - y 2 )du 2 - dx 2 - dy 2 (9.18) 



by putting 



~x = e v ' 2 x 
V = \[\V~f e~ V ' 2 x 



r = u-\{U u - V u )e M x 2 -\{U U + V u )e M y 2 

J cmi(i -ffin&u 

The wave profile in equation (9.18) now takes either of the forms 



(9.19) 



h{u) = Ci5{u) 



if ni = 2, ki = 1 



(9.20) 



9.3 The singularity structure 
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Figure 9.1 Some wave profiles for the approaching waves in the Szekeres solution 
with / = 1/2 - u n . The wavefront is unbounded if 2 < n < 4. It contains a 
step if n = 4, and is smooth if n > 6. All profiles become unbounded as u 
approaches 1. 

or 2 

HS> = 4n\^ (l - -) ^ + J ' f)3 _ yni 8(g) if m>2. (9.21) 

Particular cases in which ci = 1 and / = \ ~~ uni are illustrated in 
Figure 9.1 for ri\ = 3, 4, 5, 6 and 8. 

It may be recalled from (4.5) that = h(u). Thus, it can be clearly 
seen that the approaching gravitational wave in region II is an impulsive 
wave if ni = 2, has an unbounded wavefront if 2 < m < 4, has a step 
wavefront if n\ = 4, has a continuous wavefront if n\ > 4, and has a 
smooth wavefront if ri\ > 6. It may also be observed that, if ri\ > 4, ^ 4 
is monotonically increasing. 

9.3 The singularity structure 

In this family of solutions, there is always a scalar polynomial curvature 
singularity in region IV on the surface on which / + g = 0. This can be 
seen by computing the scalar invariant given here by 

I = 2e 2M (W 2 2 + %^° 4 ). (9.22) 



It may be noticed that an expression for the profile given by Halil (1979) 
appears to be incorrect. 
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Using (9.6) and (9.13), I can be seen to be unbounded when f + g = for 
all values of k\ and k 2 in the required range (9.9). This is similar to the 
singularity in region IV of the Khan-Penrose solution which, of course, it 
includes as a special case. 

Another singularity occurs in cases in which a wave with an un- 
bounded wave front moves in a region where the opposing wave is non- 
zero. For example, if n 2 ^ 2 and 2 < m < 4 so that 1 < k\ < |, then / is 
unbounded on the boundary u = between regions III and IV. Similarly, 
the boundary v = 0, u > is singular if rt\ > 2 and 2 < n 2 < 4 so that 
1 < k 2 < |. Apart from cases containing impulsive waves, where ri\ = 2 
or n 2 = 2, the boundaries of region IV are therefore only regular if m > 4 
for z = 1,2, which implies the further restriction | < /c| < 2 on (9.9). 
These singularities, however, are distribution-valued and their physical 
significance requires further investigation. 

The coordinate singularities in regions II and III, however, are not so 
obviously similar to the fold singularities of the Khan-Penrose solution. 
The Szekeres solutions are only flat behind the wave front in the special 
case when n\ = n 2 = 2. In the general case, regions II and III are curved 
and the analysis given in Section 8.2 cannot be applied. 

However, we may notice that in all cases with m > 2, ^4 becomes 
unbounded as / — > —1/2 in region II. This clearly indicates that, in these 
cases, the singularity in this region is a non-scalar curvature singularity. 
This has been confirmed in rigorous calculations by Konkowski and Hel- 
liwell (1989). It has also been argued in the previous chapter that these 
singularities have similar properties to the fold singularities that occur 
for colliding impulsive waves as described in Section 8.2 in terms of the 
behaviour of neighbouring geodesies. 

A similar non-scalar curvature singularity also occurs in region III if 
77-2 > 2, as ^0 then becomes unbounded as g — > —1/2. 



